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1 Introduction 

Probabilistic cellular automata (PCA) are discrete time Markov processes 
which have been intensely studied since at least Stavskaja and Pjatetskii- 
Shapiro p2] (1968). This kind of processes have as state space a product space 
X = A G where A is any finite set and G is any locally finite and connected 
graph. In this work we will focus our attention on G = Z d and A — {0, n} 
for some integer n > 1. We may regard a PCA as an interacting particle 
system where particles update its states simultaneously and independently. 
Recall that a PCA is ergodic if there exists only one invariant measure /z and 
starting from any initial measure /io the system converges to \x. 
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The aim of this paper is to use duality principles to study the ergodicity 
of two-states PCA. More precisely our work gives new sufficient ergodicity 
conditions for the expression of the PCA's local transition probabilities (see 
Theorem [5]) and show that under these conditions the invariant measure is 
shift-mixing with exponential decay of correlation. Relations between the 
PCA and the dual process (see Lemma |4] and Lemma [2] ) also allow us to 
give a very simple expression of the constant of the decay of correlation as 
a function of the radius (of the PCA) and the transition probabilities of the 
PCA (see Theorem[3|). Moreover the proof of Theorem [2] shows in detail how 
to compute the value of the invariant measure on cylinders. Results about 
the decay of correlation is an answer to a question raised in [12]. 

The existence of a dual process satisfying the duality equation (see Def- 
inition 1 and Liggett [9]) gives useful information (problems in uncountable 
sets can be reformulated as problems in countable sets) on the PCA but is 
not always sufficient to prove that a PCA is ergodic. In 12J, Lopez, Sanz 
and Sobottka introduced an extended concept of duality (see Definition 2) 
and gave general results about ergodicity (see Theorem [T]). They used this 
powerful general theory to give results on multi-state one-dimensional PCA 
of radius one and extended previous results about the Domany-Kinzel model 
(see [2] for an introduction and [7] for extensions). Previously, in [5] Konno 
has given ergodicity conditions for multi-state onc-dimcnsional PCA using 
self-duality equations. 

Even if, in some cases, the existence of null transition probabilities allows 
to prove ergodicity of a certain class of PCA (see [7] and [8]), it had been 
conjectured that in the one-dimensional case positive noise cellular automata 
are ergodic. However, P. Gacks, in 2000 , introduced a very complex coun- 
terexample (see [4] and [5]) for noisy deterministic cellular automata. In that 
case, the noisy one-dimensional cellular automata does not forget the past 
and starting from different initial distribution, the PCA may converge to 
different invariant measures. His result can be extended to noisy PCA with 
positive rates. This conjecture was formulated only in the one-dimensional 
case since in dimension 2 or higher, it is easier to show the existence of at 
least two invariant measures. For example the two dimensional Ising model 
[5] or the Toom example (see |16j ) that exhibit eroder properties. From The- 
orem [2] there exists a subclass of attractive PCA (class C) where the noisy 
conjecture is verified {p{I r ) < 1 implies ergodicity). 

In [T2], the authors explore some ergodic conditions for multi-state PCA. 
When the number of states is greater than 2, the conditions of ergodicity 
are rather restrictive in order to be able to give general results. More general 
ergodicity conditions are interesting (see [12], Section 3.2) but seems to be 
very complex when the radius grows. In this paper, we restrict the study to 
the 2 states case, which allows to show more easily general results for PCA 
of any radius. 

These sufficient ergodicity conditions can be compared to the Shlosman- 
Dobrushin condition applied to PCA (see [13] and P3]). In some examples 
(see section 3.1) our sufficient conditions induced by the concept of dual- 
ity allow to show ergodicity and decay of correlations where the Dobrushin 
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conditions can not be applied. Moreover, for some classes of ergodic PCA 
Theorem [3] gives greater constants for the decay of spatial correlation. 

This paper is organized as follows. In section 2, we present the basic 
definitions, notations and some preliminary results. In section 3, we state the 
main results, Theorem [2] and Theorem [3] We prove Theorem [2] in section 4. 
We conclude the paper in section 5 with the proof of the decay of correlation. 



2 Definitions, notations and preliminary results 

2.1 Probabilistic Cellular Automata 

We give a brief description of the theory of PCA. 

Let A be a finite set and G a countable set. Let X = A G be endowed 
with the product topology. A probabilistic cellular automata is a discrete 
time Markov process on the state space X. 

Let M(X) be the set of probability measures on X and F(X) the set of 
real functions on X which depend only on a finite number of coordinates of 



The evolution of probabilistic cellular automata is given through their 
transition operator. 

Definition 1 A local transition operator is a linear operator 



such that Pf > for all f > and PI = I, where I : X — > X is the identity 
function. 

Definition 2 A local transition operator is called independent if 



for all tp,4> G F(X) such that the finite subsets of G on which they depend 
do not intersect. 

The independent local transition operators can be defined through the values 
p x (Vi k) tor x £ G,r] £ X and k G A, as 



where for all site x G G, for all f],p x (rj,.) is a probability measure on A. 
The value p x {rj^k), called transition probabilities, represent the probability 
that the sites x G G takes the value k in the next transition if the present 
configuration of the system is n. For more details see Toom et al. J7], Maes 
and Shlosman [13] and Lopez and Sanz [TO] . 

Let d > 1 be an integer number, R a finite subset of 7L d of cardinality 
\R\ and / a map from A^ B ^ +1 to [0,1]. In the particular case G = 7/ it 



G. 



P : F(X) -> F(X) 




4 



follows from the discussion above that the discrete time Markov process 77. = 
{i]t(z) G A : t G N, z G Z d } whose evolution satisfies 



for all t G N and z G Z rf is a well defined (discrete time) stochastic process 
which from now on will be called d-dimensional PCA. Here, P stands for 
the probability measure on A z induced by the family of local transition 
probabilities. Also, let E be the expectation operator with respect to this 
probability measure. 



Let no be the initial distribution of the PCA. For any t > 0, we call n t the 
distribution of the process at time t. The measure /it is defined on cylinder 
u = N(A, </>) = {£ G A zd : £(x) = 4>{x) V x G A} for some fixed </> G A z " and 
A C Z d , \A\ < 00 by 



where Ct is the family of all cylinders of X on the coordinates of A ( the 
finite subset of Z d used to defined u). 

In this paper the notation \A\ will represent the cardinality of A when A 
is a finite subset of % d . If U = N(S, <j>) is a cylinder set, the notation \U\ will 
represent the cardinality \S\ of the set E C 7L d . In the one dimensional case we 
adopt the following notation: For any sequence of letters U = (u , ■ ■ ■ , u n ) G 
A n+1 , the set [U] s = [u ■ ■ ■ u n ] s := {x G A z \x(s) = uq, . . . , x(s + n) = u n } 
will be called cylinder and \U\ = n + 1. 

2.2 Two-state Probabilistic Cellular Automata 

In order to simplify the notation we will focus our attention on two-state 
PCA, that is to say PCA 77. on {0, 1} Z . For any positive integer r, let us 
define 

I r := {i = (ii, . . . , i d ) G Z d : -r < i u . . . , i d < r}. 
Define a family of transition probabilities {p(J) : J C I r } by 

p(J) := V{Vt+i(z) - l\Vt(z + j) = 1 : J G J}. 

Note that any PCA with state space {0, 1} Z is completely characterized 
by a positive integer number r called the radius of the PCA and the set of 
transition probabilities {p(J) : J C I r }- 

2.3 The invariant probability measure 

Definition 3 Let T be a measure-preserving transformation of a probability 
space (X, J, fi), where is the a -algebra generated by the cylinder sets on X . 
We say that the probability measure [i is T -mixing ifVU,V(z$ 



F [n t+1 {z) = a\n t (z + i) = b u V i£R} = f (a, {b t ) teR ) , 




lim n(ur\T- n v)=n(U)n(y). 



n 
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Since the cylinder sets generate the cr-algebra J*, it follows that the measure 
fi is T-mixing when the last relation is satisfied by any pair of cylinder sets 
U and V (for more details see |18|). 

2.4 Duality 

The concept of duality is a powerful tool in the theory of interacting particle 
system. It provides relevant information about the evolution of the process 
under consideration from the study of other simpler process, the dual process. 
The reformulated problems may be more tractable than the original problems 
and some progress may be achieved. Now we give the (classical) definition of 
duality taken from [9]. 

Definition 1 Let 77. and £. be two Markov processes with state spaces X and 
Y respectively, and let H (n, () be a bounded measurable function on X xY. 
The processes 77. and C. are said to be dual to one another with respect to H 
if 

E" [H( Vu C)]=E C [H(v,Ct)} 
for all 77 <E X and ( £ Y. 

Unfortunately, it is not true that every process has a dual. Recently, Lopez 
et al |12j gave a new notion of duality which extends the previous one. More 
precisely, they gave the following definition. 

Definition 2 Given two discrete time Markov processes, % with state space 
X and Ct with state space Y and H : X x Y — > R and D : Y — ► [0, 00) 
bounded measurable functions, the process 77. and £. are said dual to one 
another with respect to (H, D) if 

E no = x [H (77!, y)] = 2%)E Co=y [H (x, Ci)] . (1.1) 

2.5 Duality and sufficient conditions for ergodicity 

In order to state our results in section 3, we need to give the spirit and some 
elements of the proof of the following Theorem, which appears in [T2] , 

Theorem 1 '121 Suppose r\t is a Markov process with state space X and £t is 
a markov chain with countable state space Y , which are dual to one another 
with respect to (H,T>). If < T>(y) < 1 for all y 6 Y, then there exist a 
stochastic process £t with state space y = YU {§} with § an extra state and 
a bounded measurable function H : X x Y — > R such that 77. and £, are dual 
to one another with respect to H . Furthermore, denoting by O the set of all 
absorbing states of if 

i) the set of linear combinations of {H(.,y) : y £ Y} is dense in C (X), the 

set of continuous maps from X to R; 
ii) D(y) < 1 for any y£Q, and D := sup aeY :D( y )<i{' I %)} < 1/ 
Hi) H{.,6) = c{9) for all 9 G O with D(6) = 1; 
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then rj, is ergodic and its unique invariant measure is determined for any 
y e Y by 



(1.2) 



eeed(6)=i 



where r is the hitting time of {0 € : D{6) = 1} U {§} for £t and /t 
linit^oo fit with 



ih{y) 



H(x,y)dfit(x). 



Sketch of the proof. First recall that r is the hitting time of the dual 
process £, entering an absorbing state 9 € 0. If there exists a dual process £ 
and a function i? that satisfies the following (classical) duality equation 



E„ 



E 



io=y 



H 



(1.3) 



it is possible to show that fi s (y) = jy^H(x : y)dfi s (x) = E^ o _ [/i(^ s )]. If P{t < 
oo} = 1, it follows that 



+ hm^ooE s - o=3/ [/x(e s )|T > s]P| 0= „{t > s} 



{it=0,T< S } 



Finally, when the set of linear combinations of the set {H(.,y)\y C Z d } 
is dense in C(X) (the set of continuous functions from X to K) the sequence 
(A*n)neN converges in the weak* topology. Also, the limit measure fx does not 
depend on the initial measure /io- 

Hence, we have seen that the key point is to prove that P{r < oo} = 1. 
One way to show this, is to introduce the new type of duality (see Equation 
ll.ip . If there exists a dual process £, with state space Y that verifies the 
new concept of duality (see Equation II. ip then we can define a standard 
dual process £. with state space Y = Y U {§} and such that the set of all 
absorbing states is — U {§} where denote the set of all the absorbing 
states of £. . Here § is an extra absorbing state and the transition probabilities 
of £. satisfy 



l = So Ui = m} 



{ B(yo)P« = ? /o{6 
{ 1 - Wo) 



yi} , if 2/o,2/i 6 Y 
, if y e Y, yx = § 

. if Vq = Vx = §• 



Taking H(x,y) — H(x,y) when y G Y and H(x,§) = we obtain that the 
dual process £. satisfies the standard duality equation 11.31 and that fi(§) = 
0. Note that since D = sup ygY: D( 1 ,)<i{£ ) (2/)} < 1) a t eacn iteration the 
probability to enter the extra absorbing state § is positive and this implies 
the following result: 

Lemma 1 Under the conditions of Theorem 1, for all integer i > 1 one has 

P(r > i) < T>\ 
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Proof. By the Markov property we have that 

P t - ~ {t > i} < D x P,~ - {r > % - 1}. 
£0=3/0 L J — £0=3/0 L J 

Then, the result follows by using the mathematical induction principle. 

□ 

Note that Lemma Q] implies that P{r < 00} = 1 which finishes the proof 
of Theorem [T] 

□ 

Before stating the main results of this paper, we introduce one more 
piece of notation: let 00 1 00 denote the all one configuration, i.e. 00 1 00 = 
(l Z d(x)) 2 , gZ d- Analougsly, ^O 00 denote the all zero configuration. 

3 Main Results and Examples 

A PC A of radius r is called attractive if for any J C I r and j G I r we 
have p{ JU {j}) > p(J). We consider here the following subclass of attractive 
PC A. For any r 6 N, let p(I r ) be the set of all subsets of I r . We say that a 
two-state PCA of radius r belongs to C if its transition probabilities satisfy 
p(J) = Ylj'cJ ^(J') wnere ^ is some map from p{I r ) — > [0, 1). 

The definition of the class C is constructive. The following Proposition gives 
sufficient conditions for an attractive PCA to belong to 6. 

Proposition 1 A two-state probabilistic cellular automaton r\. belongs to C 
if its transition probabilities satisfy the following set of inequalities: 

(a) For any i <E I r , 

p({i}) >p(0). 

(b) For any 1 < k < \I r \ — f and for any jo, ■ ■. ,jk £ I r 

p({jo,...,j fc })>(-i) fc pW-E(-i) fc+1 - n E p({io,...j»». 

™=0 {! , — ,In}C{j"a. — Jh} 

Theorem 2 Let 77. &e a two-states d- dimensional probabilistic cellular au- 
tomaton of radius r that belongs to S. If p(I r ) < 1 i/ien 77. is an ergodic 
PCA and there exists a dual process £ which satisfy equation li.il Moreover, 
for any cylinder set U we can find 6 ^)k£K an d {YQe) C 1 < )keK with 
\K\ < 00 such that 

»(U) = E«* (E p £o=y(fe)te = 016-1 ^ 0} ) ■ 

fceK \;=i / 

Remark 1 ./Vote t/iat m some cases it is possible to exchange the role of the 
two states (0 <-> 1) in order to show ergodicity using the previous results. 

Corollary 1 Under the conditions in Theorem we have that if A(0) = 
then the unique invariant measure is So, where So is the Dirac measure 
on "^O 00 . Analogously, we have that if A(0) = 1 then the unique invariant 
measure p, is Si, where Si is the Dirac measure on 00 1°° . 
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Theorem 3 Lett], be a one- dimensional probabilistic cellular automaton G 6 
of radius r with p(I r ) = : D G [0, 1). Then, the unique invariant measure fi 
is shift-mixing. Also, if Ti ^ 0, for any pair of cylinders [U]o = [uq . . .Uk]o, 
[V]o = [vq ■ ■ ■ Vfe']o and t > \U\ + \ V\ we have 

|/i([f/] na- t [y] )-^([C/]o) x/i([F] )| <exp(-axi)xJf(C/,y), 

where a is the shift on {0, 1} Z , a = l/2r xln(l/D) and K(U, V) is a constant 
depending only on U , V , D and r. 

Remark 2 This last result can be extended to d- dimensional PC A. 

3.1 Examples and comparison with known results 
3.1.1 The Domany-Kinzel model 

This is a one-dimensional PCA rj, of radius r = 1 introduced in [2] with 
transition probabilities 

P{Vt+i{z) = %(z - 1, z,z + 1) = 000 or 010} = p(0) = p{{0}) = a , 
P{Vt+i(z) = l\vt(z - 1, z,z + 1) = 100 or 110} = p({-l}) 

= p({-l,0}) = oi, 
P{ry t+1 (z) = %(z - 1, z, z+ 1) = 001 or 011} - p({l}) - p({0, 1}) = oi 

and 

P{^ +1 (z) = %(z - 1, z,z + 1) = 101 or 111} = p({-l, 1}) 

= p({-l,0,l}) = O2, 

where, for any subset V G Z, ?7(V) G {0, 1} V denote the restriction of a 
configuration r] G {0, 1} Z to the set of positions in V. 

Using Proposition Q] we obtain that rj. G 6 when p({ — 1, 1}) > p({— 1}) + 
p({l}) — p(0), which is equivalent to the condition 02 > 2ai — ao. From 
Theorem [2] the PCA 77. is ergodic if p(I r ) = p({ — 1,0, 1}) = a-i < 1. From 
Theorem [3] the unique invariant measure is shift-mixing with exponential 
decay of spatial correlation such that for any pair of cylinders [U}o and [V]o 
and for all t > \U\ + \V\ we obtain 

|M[^]onff-*ry] ) -KP]o) x MMo)| < if exp(-(l/21n(l/a 2 ))i), 

where if can be explicitly computed (see the end of the Proof of Theorem[3|). 
Using Theorem [2] we can compute, for example, the measure of the cylinder 
[01]o which is 

M[oi] ) = m([i]i) - Kinh) = A({i}) - £({Q, i}) 

00 

= ^p £o={1} {6 = 0|6-i^0} 

i=l 

00 

+ 5>&={o,i}{& = 016-1 9*0}, 
*=i 

where £. is the associated dual process. 
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3.1.2 Two-dimensional example 

Let -q be a two-state, two-dimensional PCA of radius one. In this case I\ = 
{(l,k)\ — 1 < l,k < 1} is a square of 9 sites. The transition probabilities 

{p(J)\J C h} of i], are defined by p(J) — a^iio^f —ax 2^ J \ where C l k 
are the binomial coefficients. This PCA belongs to C since for any J C I\ 
we can write A(J) = a and obtain that P{J) — j'rj -M^')- This PCA is 
a kind of generalization to dimension 2 of the Domany-Kinzel model (each 
site has the same weight) with only one parameter. The sufficient ergodicity 
condition is p(I r ) < 1 which implies that a x 2 9 < 1 (a < 2~ 9 ) and the 
constant of decay of spatial correlation is a = \ ln(l/(2 9 x a)). 

3.1.3 Comparison with Dobrushin condition 

In [3], Dobrushin gives sufficient ergodicity conditions for interacting par- 
ticule systems. Using our notation, these conditions applied to PCA can be 
translated as 7 < 1 (see [T3] and [H]), where 

7 = V sup \p(J U {j}) - p(J)\. 

In the case of the Domany-Kinsel model, which belongs to the class 6, we 
obtain 7 = sup Jc/j . \p{J U {-I}) - p(J)\ + sup JcJr \p(J U {1}) - p(J)\ = 
2(a2 — ai) since ry. £ C (02 > 2ai — a ). If a 2 < 1 (condition of Theorem [2]) 
and 2(d2 — eti) > 1 the Dobrushin suficicnt conditions can not be applied. 

For the two-dimensional example we have 7 = a ^X)fe=i & x J • ^ n this case 
7 > p(I r ) and even if 7 < 1 the constant of decay of correlation | ln(l/(p(/ r )) 
is greater than i ln(l/(7)), the constant of decay of correlation given in [T5] . 
More generally, if a PCA belongs to C the sufficient condition p(I r ) < 1 
can be rewritten as p{I r ) — Sjc/ ^ ^ an< ^ ^ ne Dobrushin sufficient 

condition can be rewritten as 7 = Y^j^q JcIr A(J) x |J| < 1. 

4 Proof of Theorem [2] and Proposition [TJ 

4.1 PCA in 6 and their Dual Process 

In [12j . the authors give sufficient ergodicity conditions for one-dimensional 
multi-state PCA of radius one using a dual process satisfying equation ll.il 
Here we will use an analogous dual process to give sufficient ergodicity condi- 
tions for two-state, d-dimensional PCA of radius r using the following duality 
equation: 

E V0 = x [H(i] 1 ,Y)] = D(Y)E eo=Y H[(x,^)], (1.4) 

where 77. is a PCA with state space {0, 1} Z . The state space of the dual 
process £. is the class of all finite subsets of Z d . As in [12] we define the 
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f 1 , if x(z) = l,Vz e Y 



function H by 

^( x ' ^0 | . < >! licrwisc. 
The rule for the evolution of the process £t is given by 

where for any nonempty set J C I r we have 

P[i3(z) = {z+j|. ? e J}1 =tt(J) 



and 



= 0] = tt( 



Then, take the function D such that D(Y") = D' y ' for any finite subset 
y C Z d , where D 6 [0,1]. Note that D(0) = 1 and is the unique absorbing 
state for this dual process. 



4.2 The functions H and ft 

Note that, for this particular choice of H, we have 

ft{1 d )= [ H(x 1 Z d )d^{x)^^{ oc l co ) = 
Jx 

and 

A(0)= / #Or,0)^(*) = M{o,i} z£i ) = i, 

where X = {0, 1} Z and ^l 00 is the all one configuration (^z d ( x )) x <£Z d - The 
following Lemma is used in the proof of Theorem [2] and Theorem [3] 

Lemma 2 The set of linear combinations of {H(.,y)\y £ Z d } is dense in 
C ^{0, 1} Z , r\ the set of continuous function from {0, 1} Z to R. For any 

cylinder U = N(A,ip) C {0, 1} Z<1 (wii/i /l C Z d , \A\ < oo and ip £ A z J we 
/lave 

Y(i) 

where on £ Z, y(i) C Z d and max{|y(i)|} < oo. 

Proof. For the sake of simplicity, we only give the proof for the two-state, 
one-dimensional case. The key point of the proof consists in showing that any 
cylinder [U] t := [uo...u n ]t, (ui £ {0,1} and t,n £ N) can be decomposed 
into a non-commutative sequence of subtractions and unions of intersections 
of cylinders of the type [l] t , t £ Z. We denote by T([U] t ) this decomposition. 
One way to accomplish this decomposition is to follow the following rules: 



T([l] t ) = [l] t , T([0] t ) = {0,lf d -[l] t . 
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Then, for all t, n S Z and {/ = uq . • • u n we have 



T([C/l] t ) = T([U] t ) n [ljt+n+a. 



Thus, 



T([[/0] t ) = T([U] t ) - T([U] t ) n [l] t+2+ . 



For instance, 



r([100] )=T([10]o)-T([101]o) 

= (T([l] )-T([ll]o))-(T([10]o)n[l] 3 ) 

= [i] - [H]o - (([i]o - [H]o) n [i] 2 ) 

= ([l]o-[ll]o-([l]on[l] 2 ))U[lll] . 



Then, note that l[iooo] OJ the characteristic function of the cylinder [1000]o, 
can be written as 

l[1000] o O) = 1 [l]o( x ) + 1[111]„W - 1 [l] n[l] 2 (^) - l[U]]o( x ) 



Since for any finite subset Y C Z we have ln i6 y [i]i( a; ) = H(x,Y), it follows 
that for all n e N, t e Z and U £ {0,1}" l [u]t = J2aiH(x,Y(i)). This, in 
turn, implies that the set of linear combinations of the set {H(.,Y)\Y G Z d } 
is dense in C({0, 1} Z ). We finish the proof by observing that for any cylinder 
[17] t, we have 



Remark 3 Using the definition of H taken in \12}/ which takes into consid- 
eration the multi-state case, it is possible to prove Proposition [H for more 
general d-dimensional PC A. 

4.3 Proof of Theorem [2] 

We first establish a sequence of equalities between the transition probabilities 
of the PCA (P(J)| J € I r ) and the transition probabilities of the dual process 



H(x, {0}) + H(x, {0, 1, 2}) - H(x, {0, 2}) - H(x, {0, 1}). 




/ 5^a 4 fl-(aj ) y(*))d/i(a;) 



□ 




Hence, using the independence property of r\. we get that 



T V0=x { m {z) = 1 Vz e Y} = H P V0=x { m (z) = 1}. 
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For the left hand of equation 11.41 we have 

E So=Y [H(x,^)} = P^y{x(z) = 1 ,Vz e 
For any x £ {0, l} z<i we denote by C x the set {z £ Z d |a;(z) = 1}. Then 

w So =y{x(z) = l V* e 6} = P&=H& c Cx}- 

Using the independence property of the dual process we can assert that 

Finally we can rewrite equation 11.41 as 

= H D x P{£(z) C C x } 



z£Y zeY 

n 

z£Y 



(1.5) 



which implies that 

n = n d x a *> 

zeY zeY 

where J z = {i — z\i £ {C x D {j + z}\j £ I r }} and Z\ z is given by 



(1.6) 



A z = 7T(0) + ]T li(x(2! + *)) X ^({z}) 



X II 1 { 1 }(- T ( Z + ^))) x*"({*i> •••>**}) 
h,...,i k ei r \i=i I 



\\ l {1} (x(z + i)) \ xn(I r ). 



By simplicity of notation we write 7r(ii, ■ ■ - ik) and ■ ■ -ik) instead of 
7r({ii, . . . i k }) and p({«i, . . . ik})- 

Since equation 11.61 is true for all x £ {0,1} Z we obtain the following 
equations for 7r(.), 



p(0) = Dtt( 



p(t) = D 
P(i,j) =D 
p(i,j,k) = D 



7r(0) + 7T(t)] 

W0)+7r i) + 7r(j)+7r(i,j)] 

7r(0) + 7r(i) + 7r(j) + + k) + n(j, k) + n(i,j, k)] 



where %, j, k £ I r . 

More generally, for any < k < \I r \ — 1, 

k k-1 

p(io,...,i k ) =T>[n(9)+J2< l )+---+J2 n(lo,...,h)+ir(lo,h,...,lkj\- 

1=0 i=0l ,...,l i e{i ,—,ik} 

(1.7) 
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Since 

l-Trl 

^W + EI E *(l0,h,---h) ] +7T(J r ) = l, 

fe=o \Jo,!i,— i*e/ r 
we get that D = p{I r )- 



By definition, the dual process is completely determined by the parame- 
ters 0<7r(J)<l(Jc/ r ). From the sequence of equations 1 1.71 the dual pro- 
cess associated with the particular functions H and D exists if the transition 
probabilities of the PCA satisfy p(J) = D^/c/ 7r (^) with < p(I r ) < I r - 
In this case we have that A( J) = T)ir(J) and we claim that a PCA ?y. admits a 
dual process that satisfies the duality equation 11.11 with particular functions 
H and D given in section 4.1 if and only if this PCA belongs to the class C. 

To show that the PCA is ergodic we need to verify the three conditions of 
Theorem [TJ 

Condition i) is verified since from Lemma [2 the set of linear combinations 
of functions belonging to {#(., Y), \Y £ Z d } is dense in C ({0, 1} Z< \ m) . 
Condition ii) is satisfied since sup F _^0{D(y)} = D = p(I r ) < 1. 
Condition Hi) follows from the fact that H(., 0) = 1 and D(0) = Dl l = 1. 

Since is the only absorbing state for using Theorem Q] (equation 1.2) we 
get that for any nonemptyset Y C Z d 

A(y) = A(0)iW£t = 0} = £p&=vU* = 016-1 ^ 0}- 

t=i 

From Lemma [2l for any cylinder set U there exist S K and Y(k) finite 
subset of Jj d such that = a kjj-{Y (k)) 1 which implies the last statement 
of Theorem El 

□ 



4.3.1 Proof of Corollary^ 

When A(0) = 1, starting from any initial measure //o, we obtain that /ii = Si. 
When A(0) = 0, Theorem [2] and Lemma [2] imply that for each cylinder U 
that does not contain the point 00 0°° we get 

A*(tO = X;«i p^wft 3 0|n-i^0} J =0 
\fc=o / 

since tt(0) = 4^1 = 0. Finally we get that ^(^O 00 ) = l-/x({0, l} z<, -°°0°°) = 
1 which finishes the proof. □ 
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4.4 Proof of Proposition [T] 

Since the {tt(J)\J C I r } represent the transition probabilities of the dual 
process for all J G I r one has tt(J) > and Proposition [1] is a simple 
consequence of the following Lemma . 

Lemma 3 The transition probabilities 7r() of the dual process satisfy 

T(0) = "-f 

j) = 5 b(*i i) + - - 

■*(i,j,k) = nb(i,j,fc) -p(0) +p(i) +p(j) +p(fc) -p(i,j) -p(i,k) -p(j,k)] 
n(i,j,k,l) = n [p(M,M +p(0) -E !oe «j, M }P('o) +E { i„,, 1 }c{, J ,t,i} |, ( i|) ' il ) 
-E {io ,i 1 , i2}C { 1 , 3 ,fc,!}^°' Zl ' Z2 )l 

More generally, for any < k < \I r \ — 1 otic? /or any jo, - ■ ■ ,jk £ 



(-1)^(0) +£(-1)^' ]T pCIo,...,!,-) 

3=0 {Jo,-,ti}C{j D i-Jfc} 



Proof, of Lemma [3] From the proof of Theorem [5] a PC A belongs to class 
C if and only if the transitions probabilities pQ and 7r() satisfy the sequence 
of equations 11.71 We use mathematical induction to solve the sequence of 
equations 11.71 For the two first iterations it is easily seen that tt(0) = ^jf^, 

n(i) = p( ' ) I3 p(0) and n(i,j) = ^[p(i,j) + p(0) - p(i) ~ p{j)}- Then suppose 
that the order k is true: 



1 

D 



(-i) fe +V(0)+E(- 1 ) fc " J E pGo,...z,) 

j'=o (io,-,i,)e{jo,-,jfe} 



Using equation 11.71 we obtain that tt(Jo, ■ ■ ■ > jfe+i) equals 

k 



1 
D 



(jo,...,J fc +i)-^(0)-DE 



E 



7T(7 , . . . , ij) 



j=o \ (i ,...,i ] )e{]o,---,jk+i} 



(1-8) 

Then we suppose the rank /c true and use equation 11.81 to obtain that the 
term in p(0) in w(j , . . .,j k +i) is 



K0)-E 



E 



i=0 \io,...,ite{jo,---,jfe+i} 

K0) (-i-E c m(-i) m ) 

p(0)(- 1 + C fc+2 (-l)° + ClX 2 2 {-l) k+2 - (1 - l) fe + 2 ) 
(-i) fe+2 M0), 
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where the constants represent the binomial coefficients. Next we obtain 
that the term in J2i a e{j ,...,j k+1 } P( l o) in ""O'o, ■ ■ -,jk+i) is equal to 

-E E [ E p(Ao))(-i)' 

i=o (n ,-.i«)e{3'o,...,j*+i} \hoe{io,-,i»} / 

= e Wxx +i (-ir) 

ioe{jo,...j fc+ i} \t=o / 

= e ^o)((i-i) fc+l -^'+ i (-i) fe + i ) 

(oG{jo,--- Jfc + i} 

= E K/o)(-i) fe+1 . 

(oe{j'o,--- jfc+i} 

Note that C z fc+1 represents the number of ways to choose l\, . . . ,k in ji, . . . , jk+i 
when we have chosen lo and jo- More generally, for < M < k, the term in 
E PQo, hi) in n(j Q , . . .,j k +i) is equal to 

(io,...,(M)6{jo,--- Jfc + l} 



E E E p(h ,...,h M ) )(-ir M 

i=M (lo,-~,h)£{ja,.:,jk+i} \(ho,...JiM)£{lo,---,lj} J 

(k-M \ 
E^ +1 - M (-i)M 

E P (h,...,lM)((l-l) k+1 - M -(-l) k+1 - M ) 

(i ,...,Im )e{j'o,." 

E p ((o,..,;m)(-i) h1 - m . 

(!o,...,iA/)e{jo,--- Jfc + l} 

□ 



5 Decay of Correlation 

For the sake of simplicity we study the decay of correlation for PCA with 
state space {0, 1} Z . An extension of this result to the multi-dimensional case 
is straightforward but requires too much notation. 



5.1 Proof of Theorem [3] 

The proof of Theorem [3] requires the following two results. The second one is 
new and is a key point for the proof of Theorem [3] The first one seems to be 
well known. However, its proof can not be found or at least it is quite hard 
to be found so we provide a proof of that result. 
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Recall that \i stands for the unique invariant measure of an ergodic PCA. 

Proposition 2 Every invariant measure of an ergodic PCA is shift-invariant. 

Lemma 4 Let [U]o and [V]o be two cylinders. If /j,([U]o) — mG^Io) 
= J2 f3 l fi(B i ) and t > \U\ + \V\, then 

ii(\u] n <7-*Mo) = KP]o n [V] t ) = £ a^X*, *) Yl &A(£i)> 

where 

J2 aifj,(Ai)(*, t) J2 AA(^) := J2 a ^jK A i u i B i + *})■ 



Proof of Theorem [3] 

If D = 0, then p(0) = 0. From Corollary [1] // = 5o and has exponential 
decay of correlation. For the remainder of this proof we therefore take < 
D=p(J r )<l. 

For any finite subset E of Z and s E Z, define £J+s:={a; + s:a;€ -E}. 
We claim that for any finite subsets E and F, if t > 2Nr + \E\ + \F\ we have 

\H(E U{F + t}) fi(E) x (1(F)\ < D^ 1 ^- 

The proof of this claim uses Theorem 1 and [2] which together say that for 
any finite subset E C Z, fi.(E) = P Vo= e{i]t = 0}- This, in turn, implies that 

oo 

fi(E) =]TP, )0=£ {r = fc}, 

k=0 

where r is the hitting time for the process r\.. In fact, by Lemma [U for any 
integer N > we have 

N 

fl(E)-J2^r,o=E{r = k} 
Note that if s > 2ri + + \F\, where i is any positive integer, then 

i 

Pno=Eu{F+s}{T = i} = P V0 =e{t = i}x ^] P, ;o={F+s} {r = j} 

j=o 

i 

+ F,, 0={f+s} {t = i} x ^2f vo=Eu{f+s} {t = j}. 
It follows that if s > \E\ + \F\ + 2N x r, then 



< D 



AT+l 



1 D 
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JV 



^2^E =Eu{F+ s }{r = i} = J2 P e«=e{t = i}x ^e ={f+s}{t = i}. 



i=0 



This easily implies 



\jl{E U {F + s}) - fi{E) x fi(F)\ < D 



N+l 



1-T> 



(1.9) 



for s > \E\ + \F\ + 2N x r, which proves our claim. 

By Lemma [21 for any pair of cylinders [U]q and [V]o, there exist finite 
sequences of sets (^4^) and and finite sequences of real numbers on and 
Pi such that 



and 



Thus, by inequality 1 1.91 



aiPjfiiAi U {B. L + s} - Oip,{Ai) x ^{Bj)) < \oti0j\D 



I D 



for any pair of subsets Ai and Bj of Z and for any s > \U\ + \ V\ + 2Nr. 
It follows from this that 

I OLifa^Ai U {Bi + s}) - onfiiAi) x PjftBj) < F(U, V)T> N , 



where F{U,V) = E^kft I 1=5 • 

Using Lemma |H if t > \U\ + \ V\ we obtain 



M([C/]ona-*[y]o) -n([U]o) *tiW]o)\ <K{U,V)exp(-t 



ln(l/D) 
2~r 



where K{U, V) = F{U, V)D~(H^). 

Finally, it follows from Proposition [5] that the invariant measure is shift- 
invariant and that the exponential decay of correlations of cylinders implies 
the mixing property. 

□ 
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5.1.1 Proof of Proposition^ 

It is sufficient to show that for any cylinder [U]t, where U G {0, 1}' for some 
! 6 N , we have 

lx{a- x [U] t ) = »{[U] t ). 

Since /i is the invariant measure of an ergodic PCA rj., there exits a sequence 
(pi)i£N which converges in the weak* topology to fj,, where fii is the distribu- 
tion of a PCA rj. at time i starting from an initial distribution (iq. It follows 
that for any cylinder [U]t we have 

lim fi n {[U] t ) =K[U]t)- 

n— 'CO 

Since for any positive integer i we have 

»i{[U]t) = E Mo([^]i- i r)P noeK .] t _ ir {% 6 [U] t }, 

y 3 e{o,i}"+ 1+2i '' 

we can choose /io as a shift-invariant probability measure. Hence, for any 
positive integer i and any cylinder [U]t we have 

Hi{[U] t ) = lJii{a- l [U] t ), 

which finishes the proof. 

□ 

5.1.2 Proof of Lemma^ 

We prove the lemma using the principle of mathematical induction. First 
we prove it for the case |f7| = 1 and |V| = 1. Note that for any finite set 
{jo, ■ ■ ■ , 3k}, with j , . . . ,j k G Z, 




. . .,j k },x)dfi(x) 



= M({n[i],-|je{j ,-..,jfe}}) 

and observe that for any k > 1 e N we have 

K[i]on[i] k )=£({o}u{fc}). 

Since M ([0] fc ) = 1 - A({fc}) = A(0) - A(W) and M ([l] n [0] fe ) = fx([l] ) - 
/i([l]o n [l] fc ) = A({0}) - A({0, k}) we get, again, that 

M ([i]o n [o]*) - A({0} u 0) - A({0} u {*}). 

Furthermore, we have fi([0] o n [1]^) = A({1}) _ A({0, &})• 
Finally, note that 

M[0] o n [o] fc ) = i - m([i]o n [i] fc ) - ^([i] n [oy - /*([o] n [l]*) 

= A(0) - A({o} u {fc}) - A({0}) + A({0} u M) - A({fc}) 
+ A({0}u{fc}) 

= A(0) + A({o}u{fc})-A({o})-A(W) 
= [A(0)-A({o»] M) [A(0)-A(W)j, 
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which finishes the proof in the case \U\ = \ V\ = 1. 

Now, suppose that /z([Z7] n o- _ *[V]o) = E;,j U {S, + t}) is 

true for |t/| = |V| = n. Consider t > 2n + 1 and let [f7] , [V]o be two 
cylinders such that /x([J7]o) = E a «£(^i) an( ^ m([^]o) = E PjP-(Bj)- Since 
M[^!]o) = E <*i£(^i U {|t/|}), we get 

n([Ul]o[Vl] t ) = J2 OifrMAi U {Bj +t}U {\U\, \V\ + t}). 

hi 

Noting that /z([Vl] t ) = E a i/K{-^i + 1} U {| + 1}) we get the desired result 
for the case [Ul] H [Vl] t . 

The result for the case [Ul]o D [V0] t follows by noting that 

n(\ui] n [vo] t ) = n(pi\ n [v] t ) - n([ui] n [vi] t ) 

= £e*ft£(^U{|I/|}U{i? J -+t}) 

- J2 KSiKM u {Bj + 1] u {|c/|, |v| + 1}) 

= £ u {|t/|})(*, t) J] «iA({Bi + 0) 
-£ + t}u{|y|+*}). 

It can also be shown that 

H{[U0] n [Vl] t ) = [53a i A(^)-53a i A(i4 i U{|^|})] 
(*,t) Ea^Bi+Ou^l +*})]■ 

Finally, using that 

/x([[/o] n [K0] t ) = m([[/] n [V] t ) - n([ui] n [vo] t ) 
- fi([uo} n[vi} t ) - f,({ui} n[vi} t ) 

we can show that 

n(po] n [w] t ) = OipjiHAi u + t» 

- J2 U {Bi + t}li {\U\})}(*, t) 

- ]T ai^Ai U {-Bi + i}) U {\V\ + i}] 



which finishes the proof. 

□ 
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Final Questions 

i) Is there exist an crgodic PCA such that the unique invariant measure is 
not shift-mixing? 

ii) Is there exist an ergodic PCA such that the invariant measure has non- 
exponential decay of correlation? 
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